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INTRODUCTION 

This study describes a hybrid analog-digital pro- 
gram which, by utilizing digital timing and logical 
control signals, automatically marks time or fre- 
quency intervals on parametric plots, such as 
phase-plane or Nyquist diagrams, generated by an 
analog computer. The program is implemented on 
the TR 48 General Purpose Analog Computer with 
the DES 30 Digital Expansion System. The DES 30 
is a low-cost, general-purpose digital logic package 
which, although capable of autonomous operation for 
digital logic design or instruction, is intended to 
operate primarily with desk-top analog computers. 
With an analog computer, such as TR 48, the DES 
30 provides basic hybrid capabilities to the small 
computer facility. 

In parametric plots, the variables which are applied 
to the X and Y axes of the plotter are both functions 
of a third variable, the parameter. The parameter 
does not appear explicitly on the resulting plots, 
so that often their interpretation is complicated. 
The object of the program is to alleviate this diffi- 
culty by automatically marking segments, which 
correspond to known intervals of the parameter, 
along the solution curves (trajectories). 

Phase-plane diagrams are typical of type of plots 
being considered. Phase-plane diagrams display 
trajectories of dx/dt vs.x, where x maybe the dis- 
placement and dx/dt the velocity of a moving body. 
Both the displacement andthe velocity are functions 
of time t, the parameter. The velocity of a point 
along the phase- plane t rajectory is the phase 
velocity, ds/dt = Vx^ + y2 (not to be confused with 
the velocity of the body , dx/dt). Points at which 
ds/dt = 0 are termed singular points; at each such 
point dx/dt =dy/dt = 0 . The analysis of such phase- 
plane plots is simplified when marks, signifying 
equal time intervals, are placed on the trajectories. 



THE PROGRAM 

A simple scheme for generating such marks has 
been proposed(^). "Correction terms'* are peri- 
odically superimposed on the x and y plotter inputs 
such that a short mark, perpendicular to the slope 
of y(x), is recorded. Between marks the plotter is 
permitted to proceed without interference. In 
mathematical notation, the plotter inputs are 

X axis: x(t) — 6k ^ 
dt 

/^x M dx 
y axis: y(t) + 5k — 

where k is an arbitrary proportionality factor, and 
6 is an on-off (digital) function such that 5 = 1 for a 
short period R beginning at t=r, 2r, 3T,...nr..., and 
6=0 otherwise, r is the basic parameter interval; 
it is a constant for phase-plane diagrams, but may 
be variable for Nyquist plots. R is to be made as 
short as the dynamics of the available plotter will 
permit. The overall effect of the correction 
terms", which are added every r seconds to the 
plotter inputs, maintained for a period R, and then 
removed, is shown in Figure 1. 
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Figure 1. The Geometry of the Time Marks 
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The marks generated in this way are always per- 
pendicular to the slope of the y(x) curve, and are 
proportio nal in m agnitude to the phase velocity 
ds/dt = Vx2 + y2. Thus, on phase-plane plots, the 
magnitude of the marks diminishes as the phase 
trajectories approach a singular point; this is a 
very desirable feature since near such points 
phase contours tend to become crowded. In addi- 
tion, because of certain inherent properties of 
phase-plane contours, (dx/dt > 0 in quadrants I 
and II, limitations on the possible directions of 
contours, etc.), these marks are always on one 
side of the contour. There is no such guarantee 
for Nyquist plots. 

In the case of Nyquist plots, the derivatives dx/dt 
and dy/dt must, in general, be generated. In the 
implementation discussed here, these derivatives 
are formed digitally, although this could be done 
equally well with well-known analog programs. In 
the case of phase-plane diagrams, dx/dt = y must 
be available, and dy/dt generally is, so the situation 
is much simpler. If dy/dt is not available, either 
analog or digital differentiation programs can be 
used. 

IMPLEMENTATION WHEN DERIVATIVES MUST 
BE GENERATED 

The time derivatives are computed approximately 
as follows. The time axis is divided into equal 
intervals T. At times t = T, 2T, 3T,...nT,... the 
value of the function whose derivative is desired 
(say x(t) to be concrete) is sampled and stored 
on an analog memory (track/ store) unit. 

The output of this track/ store unit is termed 
-x(nT). At t =(n+l)T, -x(nT) is transferred to 
a second (cascaded) track/ store unit, where it 
becomes +x(nT), and a new value, -x([n+l]T), 
is sampled and stored on the first track/store. 
Forming the difference of these two track/ store 
outputs yields an increment in x: 

-Ax =x(nT) - x( [n+l] T) 

and when this increment is divided by the constant 
time increment At = T, an approximation -Ax/At 
to the derivative dx/dt is obtained. 

Figure 2a shows a possible program. Monostable 
Timer 8 (MT8) generates the signal 8 every 1 
second (r= 1 sec). The period R is controlled 
by the setting of the ON time of MT8; typically, 
R = 100 milliseconds. When § = 1, a pair of 
digitally-controlled analog switches (DA switches 
E, F) conduct, thus adding the **correction terms'' 
to the siunming amplifiers whose outputs drive the 
arm and pen of the plotter. 



Down counter 0 (DC 0) generates a carry-out (Cq) 
blip (a level which is high for one clock period 
only) every At seconds. As shown in Figure 2a, 
At = 20 milliseconds; however, in this program. 
At is adjustable from 1 to 100 milliseconds, since 
it is controlled by the setting of the preset- switch 
of DC 0, so that coarser or finer approximations 
of the derivatives can be computed (although finer 
At periods also result in increased noise sensiti- 
vity). 

The carry-out blip of DC 0 triggers MT6, which 
generates a lOO^us pulse. At the end of this lOO^us 
pulse, MT6 triggers itself again via differentiator 
7 and AND gate 13C. 

The first lOO/zs pulse generated by MT6 is applied 
to track/ stores B and D, causing them to TRACK 
for lOO^xs and return to STORE. The result of this 
operation is that the quantities -kx(nT)/At and 
+ky(nT)/At, which were stored on track/ stores A 
and C, are transferred to track/ stores B and D, 
where these quantities become +kx(nT)/At and 
-ky(nT)/At, respectively. 

The second 100 /jls pulse of MT6 is applied to 
track/ stores A and C, causing them to sample and 
store new values, -kx( [n+ 1] T)/At and ■i-ky( [n + 1] T) 
/At, respectively. At this time the summing junc- 
tion (GJ) of DA switches E and F carry the signals 
-k(Ax/At) and +k(Ay/At) respectively. These are 
the proper correction terms" except for the ob- 
vious sign inversion. 

Figure 2b is the timing diagram for the program of 
Figure 2a. This diagram shows, incidentally, that 
for all practical purposes, AND gate 13B may be 
replaced by FF 23C. 

No provision is made in this program for ad- 
justing T ; this can easily be achieved by using 
counters to frequency-divide the 1 KC source. 
For Nyquist plots with variable r (logarithmic 
frequency increments, for example), a variable- 
frequency pulse-generator may be used (2) . 

The function of the AND gate shown in broken 
lines in Figure 2a is to prevent DA switches 
E, F from conducting while the formation of the 
derivatives is in process. It may be removed at 
the cost of a '^damaged" mark now and then. 

Some saving in equipment can be effected with the 
alternative program shown in Figure 3. Here the 
function increment, Af, for the approximation to 



Reference (2) DES 20 Reference Handbook, pub. W800 2042 OA, 
EAI 1965, p. 42. 
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Figure 2a. Time Mark Program with Derivative Computation 
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Figure 2b. Timing Diagram 
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the derivative, is formed by subtracting a sampled 
and stored past value of the fimction from the 
present value, but no attempt is made to store 
the present value. Thus the **present" value 
actually changes during the time that the mark 
is recorded, resulting in marks which are not 
so well-defined but are nevertheless quite ac- 
ceptable for many applications. This program dif- 
fers from the one in Figure 2a also in that deri- 
vatives are formed only prior to the generation 
of the mark, rather than continually. Here At is 
controlled by MT7 while MT8 generates the signal 
d (controls R). 

IMPLEMENTATION WHEN DERIVATIVES ARE 
AVAILABLE 

In the case of phase plane plots, dx/dt =y is avail- 
able and so is, generally, dy/dt. Hence the program 
of Figure 4 suffices. Again, control over r can be es- 
tablished with counters. The MT shown controls R. 



PROGRAMMING INFORMATION 

Attenuators which feed digital devices such as 
track/ stores or DA switches must be set with 
the device in the ON (conducting, TRACK) state, 
since these devices present different loads in the 
ON and OFF states. 

The program is initialized by setting the proper 
ON time controls of monostables and preset 
switches of down counters. Note that At is set on 
several attenuators, in addition to a monostable. 
Then clear the DES 30, place it in 1 MC RUN, and 
put the TR 48 from IC to OP. 

The program of Figure 2a requires 3 Electronic 
Comparator 40.488 units on the TR 48, and 3 MT- 
DIF, 1 GPR, 1 AG trays on the DES 30. The pro- 
gram of Figure 3 requires 2 40.488 Comparator 
units on the TR 48 and the same complement of 
equipment on the DES 30 as in Figure 2a. 
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Figure 3. An Alternative Time Mark Program 
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Figure 4. Time Marks for Phase Plane Diagrams 
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